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ABSTRACT 

X-ray observations of blazars associated with the OVV (Optically Violently 
Variable) quasars put strong constraints on the e + e~ pair content of radio-loud 
quasar jets. From those observations, we infer that jets in quasars contain 
many more e + e~ pairs than protons, but dynamically are still dominated by 
protons. In particular, we show that pure e + e~ jet models can be excluded, 
as they overpredict soft X-ray radiation; likewise, pure proton-electron jets 
can be excluded, as they predict too weak nonthermal X-ray radiation. An 
intermediate case is viable. We demonstrate that jets which are initially 
proton-electron ("proto-jets") can be pair-loaded via interaction with 100 - 300 
keV photons produced in hot accretion disc coronae, likely to exist in active 
galactic nuclei in general. If the coronal radiation is powered by magnetic flares, 
the pair loading is expected to be non-uniform and non-axisymmetric. Together 
with radiation drag, this leads to velocity and density perturbations in a jet 
and formation of shocks, where the pairs are accelerated. Such a scenario can 
explain rapid (time scale of ~ a day) variability observed in OVV quasars. 



Subject headings: galaxies: jets — plasmas — radiation mechanisms: 
non-thermal 
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1. INTRODUCTION 

One of the basic unresolved questions regarding the nature of jets in radio loud quasars 
is that of their composition: are they made from protons and electrons, or electron-positron 
pairs, or from a mixture of both? Arguments in favor of proton-electron jets in quasars 
have been recently advanced by Celotti & Fabian (1993). Using synchrotron self-Compton 
constraints from radio-core observations and information about energetics of jets from 
radio-lobe studies, those authors showed that in the case of pure electron-positron jets the 
required number of e + e~ pairs is too high to be delivered from the central engine. The limit 
is imposed by the annihilation process (Guilbert, Fabian, & Rees 1983). 

On the other hand, the recently discovered circular polarization in the radio cores 
of the 7-ray bright OVV quasar 3C 279 and several other objects and its interpretation 
via the "Faraday conversion" process suggests that the jet plasma is dominated by e + e~ 
pairs (Wardle et al. 1998; Wardle & Homan 1999). The fact that jets are likely to be pair 
dominated has also been inferred from synchrotron self-Compton analyses of compact radio 
components in radio galaxy M87 (Reynolds et al. 1996), and in quasar 3C 279 (Hirotani et 
al. 1999). In this paper we derive constraints imposed on the pair content of quasar jets by 
X-ray observations of OVV quasars, i.e. those radio loud quasars which are observed at 
very small angles to the jet axis, and often detected in the MeV - GeV 7-ray regime. Our 
results suggest that the pair content of quasar jets is high, but that dynamically the jets 
are dominated by protons. 

The question of composition of the jet plasma is closely related to that of the formation 
of the jet. Jets can be launched as outflows dominated by Poynting flux, generated in the 
force-free magnetosphere of the black hole, or as hydromagnetic winds driven centrifugally 
from an accretion disc (see, e.g., a review by Lovelace et al. 1999). Electromagnetically 
dominated outflows are converted to pair dominated jets (Romanova & Lovelace 1997; 
Levinson 1998), whilst hydromagnetic winds give rise to proton-electron dominated jets 
(Blandford & Payne 1982). While the pair-dominated jets are predicted to be relativistic, 
the proton-electron jets can be either relativistic or non-relativistic, depending on whether 
the magnetic forces dominate over gravity in the accretion disc corona (Meier et al. 1997). 
In particular, relativistic hydromagnetic jets can be launched deeply in the ergosphere 
of fast-rotating black holes (Koide et al. 1999). If this is the case, and the surrounding 
accretion disc corona is as hot as inferred from the spectra of Seyfert galaxies (see, e.g., 
Zdziarski et al. 1994; Matt 1999), the proton-electron jets can be pair- loaded via the 
interactions with the coronal hard X-rays / soft 7-rays. We demonstrate that this process 
is efficient enough to provide the number of pairs required to account for the observed X-ray 
spectra of OVV quasars. Furthermore, the rapid X-ray variability of Seyferts (Green et al. 
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1993; Hayashida et al. 1998) indicates that the corona is likely to have dynamical character 
(as would be if it is powered by magnetic flares), and thus the hydromagnetic outflows are 
expected to be loaded by pairs nonuniformly and non-axisymmetrically. This, together with 
the radiation drag imposed by the coronal and disc radiation fields on pairs, can lead to a 
modulation of the velocity and density of the plasma in a jet, and therefore to production 
of shocks and acceleration of particles responsible for the variable nonthermal radiation. 

Our paper is organized as follows: In §2 we demonstrate that pure pair jets overpredict 
the soft X-ray flux; in §3 we show that jets which are dynamically dominated by protons 
must be heavily pair-loaded in order to provide the observed nonthermal X-ray radiation; 
and in §4 we show that hydromagnetic outflows can be pair-loaded due to interaction with 
hard X-rays / soft 7-rays from the hot accretion disc corona. Our results are summarized 
in §5. 



2. ELECTRON-POSITRON JETS 



2.1. Cold pairs 



If the relativistic jet is dynamically dominated by cold e + e~ pairs, then the external 
UV photons will be Comptonized by those pairs and thus boosted in frequency by a square 
of a bulk Lorentz factor T and collimated along the jet axis. In this case, in addition to 
the nonthermal radiation from the jet - which results in the phenomenon called blazar - 
the observer located at o b s < 1/r should see a soft X-ray "bump" superimposed on the 
continuum in OVV quasars. Such spectral feature, predicted theoretically by Begelman & 
Sikora (1987), has not been observed, and this fact can be used to derive an upper limit for 
a pair content of quasar jets (Sikora et al. 1997). 



Luminosity of the soft X-ray bump produced by the above "bulk-Compton' : 
process is 
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A is the beaming amplification factor, Ld is the luminosity of an accretion disc, £ is the 
fraction of the accretion disk which at the given distance r is isotropized due to rescattering 
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or reprocessing, n e is the density of electrons plus positrons, dV = 7ia 2 dr, Tj = n e aaT, a 
is the cross-section radius of a jet, and r is the distance at which the jet is fully formed 
(accelerated and collimated). 

Assuming that jet is conical and that at r > r$ the electron/positron number flux is 
conserved (no pair production), we have n e oc 1/r 2 , and Tj = ri/r, where by r\ we denote 
the distance at which Tj = 1. Then, for jets with a half-angle 9j = a/r < 1/r, Eqs. (|l])-([f) 
give 

1, w x , fln a + l if r 7 (r ) > 1 
3 V JAS 7 I £ if Tj r < 1 ' 



J BC 



where we used A = T 2 . The value of ri depends on the electron/positron number flux, 
dN e /dt, which - for energy flux in a jet Lj dominated by kinetic energy flux of cold pairs - 
is equal to Lj/m e c 2 T. Noting that dN e /dt ~ n e ciia 2 , we obtain 

1 a TLj 17 Lj Ag , . 

r% = — ~ ^- — ~ 10 — p — cm , (5) 

o^n e Qj 'am e c J '9jT 8jT 

If ro > ri then n should be treated only as a formal parameter which provides normalization 
of Tj. However, noting the very large value of r\ one can expect that r < r\ and, therefore, 
the predicted Bulk-Compton luminosity is 

L BC > 3 x 10 47 (r^)(eL d ) 4 5(r/10) 3 (ln(r 1 /r ) + l^rgs" 1 . (6) 

The BC spectral component peaks at hv ~ Y 2 hvuy ~ (T/10) 2 keV, where typical 
luminosities observed in OVV quasars are ~ 10 46 erg s" 1 (Sambruna 1997), while the 
spectra are consistent with simple power laws (cf. Kubo et al. 1998). Thus, one can 
conclude that pure electron-positron jet models can be excluded as overpredicting soft 
X-ray radiation of OVV quasars. 



2.2. Relativistic "thermal" pairs 

It is of course possible that because of inefficient cooling of electrons/positrons 
below a given energy, the multiple reacceleration process balances adiabatic losses in the 
conically diverging jet, and the pairs, once accelerated, remain relativistic forever. If the 
relativistic electrons are narrowly distributed around some 7, then taking into account 
that \dE e /dt\ oc T 2 7 2 , n e oc Lj/j, and 7*1 ~ 77(7 = l)/7, and assuming r > 7*1, the 
bulk-Compton luminosity would be 

L BC = 3 x 10 47 (^,)(^) 45 (r/10) 3 (r 1 /r )7 * 3 x 10 48 ^M^Mlh^ ergs -i , (7) 



- 5 - 



and would peak at ~ his ~ T 2 ^ 2 u uv ~ 7 2 (T/10) 2 keV. No such "bumps" have been detected 
at keV energies. Alternatively, one can speculate that the X-ray spectra of OVV quasars 
consist of superposed multiple "thermal" peaks produced over several decades of distance 
(cf. Sikora et al. 1997). In this case it may be possible to match the observed X-ray 
spectral slopes, but nonetheless, the model predicts too large luminosity. 



The other extreme would have no e + e~ pairs in a jet. For a given energy flux in the 
jet, Lj, which now is proportional to n p m p , the number of electrons in a jet is m e /m p 
times smaller than the number of electrons plus positrons in the jet made from cold pairs. 
Thus, noting that ri(n e = n p ) = (m e /m p )ri(n p = 0) ~ 0.5 x 10 14 cm, one can find that 
the proton-electron jets do not overproduce the soft X-ray luminosities observed in OVV 
quasars, provided that r > 15 r! ~ 10 15 cm. 

However, such pure proton-electron jets are relatively inefficient in producing the 
nonthermal radiation, and this is for the same reason as above - low number of electrons. 
This is apparent from a study of the low energy tails of the nonthermal radiation 
components, where the requirement for the number of electrons is largest. In the case of 
synchrotron radiation, such tails are not observed because they are self-absorbed, and the 
only spectral band where the presence of the lower energy relativistic electrons can be 
evident are the soft and mid-energy X-rays, 0.1 - 20 keV. 

It has been shown in many papers (see, e.g., Sikora, Begelman, & Rees 1994) that 
the 7-ray spectra observed by the EGRET instrument in OVV quasars are likely to be 
produced by Comptonization of external diffuse radiation field, via the so-called external 
radiation Compton (ERC) process. Can this process be also responsible for the X-ray 
spectra of OVV quasars? In the one zone model and for the narrow spectral distribution of 
the soft radiation field, the power law X-ray spectra are produced by electrons with energy 
distribution obeying n'^ = C n , ~f~ s , where s = 2a x + 1 (note that throughout this paper, 
all quantities except for 7 are primed if measured in the frame co-moving with the jet). 
Assuming that at a distance r^, where the blazar phenomenon is produced, all available 
electrons are accelerated and that energy flux in the jet is dominated by cold protons, i.e., 
that 
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The ERC luminosity is then given by 
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and 



^(4/3)ry^t. (13) 

The external radiation field in quasars, as seen in the jet frame at a distance 

a ctfiT 17 (t f/ /ld)(I7l0) 2 , , 

r/i = j - -jj- - 2-6 x 10 17 1 fW ^ ' cm, (14) 

where tfi is the time scale of duration of flares observed in OVV quasars, is dominated by 
two components, broad emission lines and near infrared radiation from hot dust. The broad 
emission lines provide radiation field with hv ext ~ 10 eV and 

= 3 ' 9xl ° (wn«i7W ergcm ' (15) 

The spectrum of the infrared radiation from hot dust, on the other hand, peaks around 
hu ext = 3kT ~ 0.26(T/1000) eV, and at a distance rfi < r IR = ^ L d / AucjsbT 4 , has energy 
density 

u diff{IR) * tm^f-T* ~ 7.6 x lO- 3 ^ (tJ^) erg cm" 3 , (16) 

where T is the temperature of dust, <7sb is the Stefan-Boltzman constant, and is the 
fraction of the central source covered by the innermost parts of a dusty molecular torus. 
For T ~ 10 and typical 1-20 keV spectral index a x — 0.6 (Kii et al. 1992; Kubo et al. 
1998), Comptonization of broad emission lines gives 

(L v u) ciBEL) - 6.4 x 10 43 |^ {^y^L^m ergs" 1 , (17) 



while Comptonization of near infrared radiation gives 



Mcm K ,1 x 10- (|f (J^L) (*) ^ j ^„L,«er gS - . (18) 
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Note that in the case of Comptonization of UV photons, radiation at 1 keV is produced 
by electrons which are only weakly relativistic, with 7 ~ 1 (see Eq. and therefore this 
implies that 7 m j„ also needs to be ~ 1. In the case of Comptonization of near infrared 
radiation, 7 min < 50/T. 



As one can see from Eqs. (|17D and (jig) , Comptonization of external radiation 
by relativistic electrons in the proton-electron jets gives 1 keV luminosities which are 
~ 100/Lj i46 times smaller than observed. 

Let us now determine whether the observed X-ray luminosities can be produced by the 
pure proton-electron jets via the synchrotron self-Compton (SSC) process, in addition to 
the ERC process responsible for the hard 7-ray emission. Luminosity of the SSC radiation 
can be estimated using the formula 



(L u u) ssc ~ T\L v ,v') S sc ^ T\^N yl )m e c 2 
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dt' 
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and 



syn ~ 2vrca 2 r 4 ' 1 } 



v ~ (4/3) 7 2 ^n, m , (22) 

where /i^ S2/ „. m ~ 0.1 eV is the typical location of the synchrotron spectrum peak in OW 
quasars (cf. Fossati et al. 1998). 



As it is apparent from Eq. (|2"2|), production of IkeV radiation by SSC process involves 
electrons with 7 ~ 100. Therefore, 7 min is not restricted to such low values as in the case 
of the ERC processes, and, in principle, for -y m in ~ 100, the SSC model can reproduce the 
observed soft X-ray luminosities: 

^ 73 x w,s ^mm (tSv) 01 (wioo)"-g^ , (23) 

where as before we used T = 10 and ax = 0.6. However, since the electrons which 
produce X-ray spectra via the SSC process have the same energy range as those electrons 
which produce 7-rays above 1 MeV, the spectral slopes of both should be similar. The 
observations show that this is not the case; the 7-ray spectra above 1 MeV are much 
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steeper than the X-ray spectra in the 1-20 keV range, typically by Aa ~ 0.5 (Pohl et 
al. 1997). This contradiction can be eliminated by assuming that most of the electrons are 
injected with 7 > 500. The X-rays are then produced by electrons which reach energies 
appropriate for X-ray production (100 < 7 < 500 for 1-25 keV) by radiative energy losses 
and the resulting slope is «i ~ 0.5 (Ghisellini et al. 1998; Mukherjee et al. 1999). 

Summarizing this section we conclude that: 

• Production of hard X-ray spectra by ERC process requires 7 miri < 10 and the pair to 

proton number ratio 

~50 LsX ' 46 . (24) 
n p Lj^ 

where L S x ~ 10 46 erg s _1 is the typical luminosity observed in OVV quasars around 1 
keV; 

• Hard X-ray spectra can be produced by pure proton-electron jets via the SSC process, 

but this requires extremely high values of minimum electron injection energies. 



4. PAIR PRODUCTION AND VARIABILITY 

We propose a scenario where jets are launched as proton-electron outflows in the 
innermost parts of the accretion flow and are loaded by pairs due to interactions with hard 
X-rays / soft 7-rays produced in the hot accretion disc coronae. This can well occur via a 
two-step process. The first step is Compton boosting of coronal photons (with initial energy 
of 100 - 300 keV) up to few MeV by cold electrons in the outflow propagating through the 
central region (Begelman & Sikora 1987). Provided that luminosity of the coronal radiation 
at > 100 keV is L S1 ~ 10 46 erg s _1 , as can be deduced from extrapolation of 2 — 10 keV 
spectra observed in non-OVV radio- loud quasars (see, e.g., Cappi et al. 1997; Xu et al. 
1999), one can find that opacity for the above interactions is very high, 

(L^/lO^ergs- 1 ) 

r e7 ~ n x r corona a T ~ 15 ^^^—^ ~ 10 i5 cm) ■ ( 25 ) 

This means that each electron in the inner parts of the outflow (essentially forming a 
"proto-jet") produces on the order of 10 or more 1-3 MeV photons. The second step is 
the absorption of MeV photons by the coronal (100 - 300 keV) photons in the pair creation 
process. The pairs created in such a manner are dragged by the jet, but before leaving the 
compact coronal radiation field, they produce a second generation of MeV photons, and 
they in turn produce next generation of pairs. Such pair cascade can continue until the time 
when the proto-jet becomes opaque for coronal radiation, i.e., when n e r corona <jT ~ 1- Within 
this limit, the electron/positron flux integrated over the cross-section of the proto-jet can 
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reach the value 

N e ~ n e dl,ir 2 corona ~ ; (26) 

where f2j is the initial solid angle of the outflow. Comparing this electron/positron flux 
with the total proton flux 

N p ~ -M~ 2 , (27) 
1 m p c z 

we find that proton-electron winds can be loaded by pairs in the central compact X-ray 
source up to the value 

repairs TYlpC T corona VliJ^ r corona f2j T , no , 

— - — ~ — — ~ — ~ 30 — — . (28) 



n 



p 



N p 2a T Lj 3 x 10 15 cm L jA6 3 



This corresponds roughly to the pair content given by Eq. (|4j) , provided that Q $ is not very 
small. Note that a large initial opening angle of the central outflow is expected, as jets with 
Lj ~ 10 46 erg s _1 carry too much momentum to be effectively collimated by the innermost 
parts of the accretion disc corona. Due to radial quasi-expansion, the outflow is rapidly 
diluted and can be collimated to the narrow jet by disc winds at r > 100 gravitational radii 
(cf. Begelman 1995). 

It should be mentioned here that loading of quasar jets by pairs via absorption of 
7-rays produced within the jet by external radiation field has been also proposed by 
Blandford and Levinson (1995) (BL95). However, their scenario is very different from ours 
in many respects. In the BL95 model, both pairs and nonthermal radiation are produced 
over several decades of distance; in our model, pair production is taking place at the base 
of a jet, whereas nonthermal radiation is produced at distances 10 17 - 10 18 cm. The BL95 
model involves relativistic electrons, and pairs are produced by absorption of nonthermal 
radiation extending up to GeV energies, while our pair loading scenario involves cold 
electrons (as measured in the jet comoving frame), and pairs are produced by absorption 
of photons with energy 1-3 MeV. In their scenario, Comptonization of the UV bump by 
relativistic pair cascades leads to a production of a power-law X-ray spectrum which is 
softer than that observed in OVV quasars; in our scenario — in the region where pairs are 
injected, i.e., at the base of a jet — the UV bump is Comptonized only by cold pairs, and 
this leads to a production of radiation only around hv V y{F /3) 2 ~ 100 eV. Due to the wide 
opening angle of a jet at its base, this radiation is much less collimated than nonthermal 
radiation produced at larger distances, and therefore in OVV quasars, it may be relatively 
inconspicuous. The 100 eV excess can be detectable eventually in steep spectrum radio 
loud quasars, which have jets pointing further away from our line of sight. However, due to 
absorption by the ISM in the host, or our own Galaxy, this excess is predicted to be weak, 
and difficult to detect. 
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Another attractive feature of our model is that the pair loading of a jet via interactions 
of the proto-jet with the hard X-rays / soft 7-rays produced by accretion disk coronae 
can be responsible for fast (~ day) variability observed in OVV quasars. We know from 
observations that the (presumably isotropic) X-ray emission from Seyfert galaxies - and 
thus, by analogy, the non-jet, isotropic component in radio loud quasars - is rapidly variable 
(although in OVV quasars, this component is "swamped" by stronger, relativistically 
boosted flux). This suggests that the corona is likely to have dynamical character: it may 
be powered by magnetic flares, or else by a possible instability of the innermost region of 
the accretion disk. In either case, the jet is expected to be loaded by pairs non-uniformly 
and non-axisymmetrically. The patches of the local pair excesses in a jet suffer large 
radiation drag (Sikora et al. 1996) and are forced to move slower than the surrounding gas. 
Therefore, they provide natural sites for shock formation and particle acceleration. 

While the above mechanism is viable as an explanation of rapid X-ray and 7-ray 
variability observed in OVV quasars, we note that pair density variations, as modulated 
by magnetic flares in the disk, are too rapid to produce variability of the radio flux. The 
long-term (months to years) variability in OVV quasars observed in all spectral bands 
including radio, are more likely to result from modulation of the variable flux of protons. 
Such modulation can be induced by the variability of the accretion rate in the inner parts 
of the accretion disc. The observed long term optical variability in both radio-loud and 
radio-quiet quasars supports this view (see, e.g., Giveon et al. 1999 and references therein). 

5. SUMMARY 

• Models of quasar jets consisting purely of e + e~ pairs can be excluded because they predict 

much larger soft X-ray luminosities than observed in OVV quasars. On the other 
hand, models with jets consisting solely of proton-electron plasma jets are excluded, as 
they predict much weaker nonthermal X-ray radiation than observed in OVV quasars. 
Spectra of nonthermal flares in those objects can be explained in terms of a simple 
homogeneous ERC model, provided that the number of pairs per proton reaches values 
-50 (Lsx/Lj). 

• We suggest that initially, jets consist mainly of proton-electron plasma (where the protons 

provide the inertia to account for the kinetic luminosity of the jet), and subsequently are 
loaded by e + e~ pairs by interactions with hard X-rays / soft 7-rays from hot accretion 
disc coronae. This requires that the coronal temperatures reach values ~ 100 keV, 
which are consistent with observations of Seyfert galaxies, with spectra uncontaminated 
by relativistic jets. 

• Non-steady and non-axisymmetric pair loading of jets by X-rays from magnetic flares 

in the corona can be responsible for short term (~ day) variability observed in OVV 
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quasars. It should be emphasized here that alternative mechanisms of variability, such 
as modulation of the total energy flux in a jet by accretion rate or precession of a jet 
(cf. Gopal-Krishna & Wiita 1992) cannot operate on such short time scales. The lack 
of rapid (time scale of ~ day), high amplitude variability in the UV band of radio lobe 
dominated quasars supports this view. 
• The dissipative sites in quasar jets, where electrons/positrons are accelerated and 
produce nonthermal flares observed in OVV quasars, can be provided by shocks 
produced by collisions between inhomogeneities induced by non-uniform pair loading 
of the proto-jets. These shocks (and therefore particle acceleration) can be amplified 
eventually at a distance 0.1 - 1 pc due to reconfinement of a jet by the external gas 
pressure (Komissarov & Falle 1997; Nartallo et al. 1998). 

We are grateful to Annalisa Celotti, the referee, whose thoughtful comments helped 
us to substantially improve the paper. M.S. thanks Annalisa Celotti and Paolo Coppi for 
stimulating discussions, and, in particular, for pointing out that the pair-loading process 
in proto-jets can saturate. We would like to thank the Institute for Theoretical Physics at 
U.C. Santa Barbara for its hospitality. This project was partially supported by ITP/NSF 
grant PHY94-07194, NASA grants and contracts to University of Maryland and USRA, 
and the Polish KBN grant 2P03D 00415. 
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